This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

**Table 1.** Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue	Howard	Bayes
Time series	1999 - present	1977 - present
Bias in SWHS	Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests.	Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases	Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data.	Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations	Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas.	Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation	Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with.	By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Results

Model: LB_fit_3pH
Model Run: Gen4int_indcomp_swhsR_FULL_pHB4pars_re0d_thru2024_2e+05_2025-10-30 (Naming convention: modelname_thruYEAR_iters_rundate)
Proportion of parmaters converged (Rhat <= 1.1): 0.9947
See end of document for summary of unconverged parameters.

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

And black rockfish releases…

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Total Biomass Removal Estimates

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Model fit

Logbook residuals

Figure 20.- Residuals from logbook harvests.

SWHS residuals

Figure 21.- Residuals from SWHS harvests.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

Figure 23.- Mean percent of harvest by charter anglers.

When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(slope|non-pelagic & non-yellowye) For release estimates

## NULL

Weight Fits

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.2419355

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.12 while 0.88 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable	area_n	year	user	mean	median	lower	upper	Rhat	sd	skew	area
p_pelagic	11	1	1	0.474	0.468	0.237	0.728	1.026	0.127	0.047	CSEO
p_pelagic	11	1	2	0.409	0.403	0.192	0.651	1.011	0.122	0.049	CSEO
p_pelagic	11	10	1	0.477	0.476	0.233	0.724	1.022	0.126	0.008	CSEO
p_pelagic	11	10	2	0.412	0.403	0.189	0.668	1.022	0.124	0.073	CSEO
p_pelagic	11	11	1	0.471	0.468	0.242	0.713	1.021	0.122	0.025	CSEO
p_pelagic	11	11	2	0.409	0.402	0.191	0.668	1.022	0.125	0.056	CSEO
p_pelagic	11	12	1	0.475	0.472	0.237	0.717	1.015	0.127	0.024	CSEO
p_pelagic	11	12	2	0.406	0.398	0.192	0.650	1.017	0.121	0.066	CSEO
p_yellow	11	12	1	0.518	0.517	0.235	0.802	1.011	0.149	0.007	CSEO
p_pelagic	11	13	1	0.476	0.476	0.232	0.730	1.014	0.128	0.000	CSEO
p_pelagic	11	13	2	0.409	0.404	0.187	0.655	1.015	0.123	0.041	CSEO
p_pelagic	11	14	1	0.485	0.485	0.252	0.730	1.028	0.125	0.000	CSEO
p_pelagic	11	14	2	0.414	0.408	0.193	0.671	1.018	0.123	0.049	CSEO
p_pelagic	11	15	1	0.444	0.440	0.216	0.705	1.030	0.127	0.032	CSEO
p_pelagic	11	15	2	0.403	0.396	0.187	0.669	1.015	0.124	0.056	CSEO
p_pelagic	11	16	1	0.491	0.489	0.258	0.723	1.022	0.122	0.016	CSEO
p_pelagic	11	16	2	0.414	0.409	0.185	0.665	1.025	0.125	0.040	CSEO
p_pelagic	11	17	1	0.481	0.477	0.256	0.714	1.011	0.121	0.033	CSEO
p_pelagic	11	17	2	0.410	0.404	0.187	0.669	1.022	0.123	0.049	CSEO
p_pelagic	11	18	1	0.481	0.480	0.245	0.726	1.019	0.123	0.008	CSEO
p_pelagic	11	18	2	0.411	0.410	0.187	0.658	1.010	0.122	0.008	CSEO
p_pelagic	11	19	1	0.494	0.494	0.262	0.728	1.022	0.121	0.000	CSEO
p_pelagic	11	19	2	0.416	0.408	0.190	0.675	1.018	0.126	0.063	CSEO
p_pelagic	11	2	1	0.474	0.473	0.242	0.719	1.019	0.127	0.008	CSEO
p_pelagic	11	2	2	0.410	0.404	0.188	0.664	1.011	0.125	0.048	CSEO
p_yellow	11	2	2	0.426	0.413	0.169	0.737	1.011	0.150	0.086	CSEO
p_pelagic	11	20	1	0.488	0.483	0.263	0.733	1.029	0.123	0.041	CSEO
p_pelagic	11	20	2	0.416	0.411	0.193	0.672	1.017	0.124	0.040	CSEO
p_pelagic	11	21	2	0.415	0.409	0.189	0.678	1.011	0.125	0.048	CSEO
p_yellow	11	21	2	0.428	0.419	0.162	0.750	1.013	0.154	0.058	CSEO
p_pelagic	11	22	2	0.417	0.411	0.196	0.673	1.016	0.123	0.049	CSEO
p_yellow	11	22	1	0.521	0.519	0.242	0.809	1.014	0.151	0.013	CSEO
p_pelagic	11	23	2	0.452	0.449	0.233	0.696	1.012	0.121	0.025	CSEO
p_yellow	11	23	1	0.537	0.536	0.268	0.812	1.011	0.143	0.007	CSEO
p_yellow	11	25	1	0.486	0.478	0.222	0.776	1.011	0.145	0.055	CSEO
p_yellow	11	25	2	0.442	0.432	0.170	0.758	1.011	0.154	0.065	CSEO
p_yellow	11	26	1	0.512	0.509	0.243	0.794	1.011	0.145	0.021	CSEO
p_yellow	11	28	1	0.486	0.481	0.214	0.774	1.011	0.149	0.033	CSEO
p_yellow	11	29	1	0.501	0.499	0.230	0.800	1.014	0.148	0.014	CSEO
p_pelagic	11	3	1	0.475	0.471	0.237	0.724	1.021	0.128	0.031	CSEO
p_pelagic	11	3	2	0.409	0.403	0.182	0.671	1.013	0.126	0.048	CSEO
p_pelagic	11	4	1	0.479	0.478	0.248	0.723	1.013	0.125	0.008	CSEO
p_pelagic	11	4	2	0.407	0.401	0.187	0.650	1.017	0.121	0.050	CSEO
p_yellow	11	4	2	0.428	0.418	0.175	0.738	1.016	0.147	0.068	CSEO
p_pelagic	11	5	1	0.476	0.474	0.238	0.720	1.011	0.127	0.016	CSEO
p_pelagic	11	5	2	0.411	0.402	0.189	0.674	1.021	0.124	0.073	CSEO
p_pelagic	11	6	1	0.481	0.479	0.247	0.730	1.014	0.126	0.016	CSEO
p_pelagic	11	6	2	0.412	0.408	0.185	0.667	1.023	0.124	0.032	CSEO
p_yellow	11	6	1	0.516	0.516	0.238	0.806	1.011	0.151	0.000	CSEO
p_pelagic	11	7	1	0.476	0.475	0.233	0.722	1.022	0.126	0.008	CSEO
p_pelagic	11	7	2	0.411	0.407	0.192	0.666	1.023	0.123	0.033	CSEO
p_pelagic	11	8	1	0.478	0.479	0.241	0.731	1.015	0.128	-0.008	CSEO
p_pelagic	11	8	2	0.407	0.401	0.189	0.654	1.019	0.123	0.049	CSEO
p_pelagic	11	9	1	0.474	0.470	0.237	0.727	1.017	0.128	0.031	CSEO
p_pelagic	11	9	2	0.411	0.409	0.189	0.662	1.024	0.125	0.016	CSEO
p_pelagic	13	1	1	0.343	0.328	0.123	0.634	1.074	0.135	0.111	NSEI
p_pelagic	13	1	2	0.149	0.132	0.043	0.346	1.066	0.079	0.216	NSEI
p_yellow	13	1	1	0.430	0.374	0.129	0.982	1.029	0.224	0.250	NSEI
p_yellow	13	1	2	0.431	0.380	0.127	0.981	1.028	0.223	0.228	NSEI
p_pelagic	13	10	1	0.352	0.341	0.125	0.637	1.060	0.132	0.083	NSEI
p_pelagic	13	10	2	0.155	0.138	0.046	0.357	1.074	0.081	0.211	NSEI
p_yellow	13	10	1	0.429	0.381	0.125	0.982	1.032	0.224	0.214	NSEI
p_yellow	13	10	2	0.433	0.384	0.126	0.984	1.029	0.224	0.219	NSEI
p_pelagic	13	11	1	0.354	0.340	0.134	0.647	1.053	0.133	0.105	NSEI
p_pelagic	13	11	2	0.154	0.138	0.046	0.348	1.049	0.079	0.202	NSEI
p_yellow	13	11	1	0.433	0.380	0.126	0.982	1.028	0.225	0.236	NSEI
p_yellow	13	11	2	0.434	0.384	0.129	0.983	1.031	0.224	0.224	NSEI
p_pelagic	13	12	1	0.360	0.350	0.140	0.632	1.048	0.131	0.076	NSEI
p_pelagic	13	12	2	0.154	0.141	0.045	0.340	1.064	0.078	0.167	NSEI
p_yellow	13	12	1	0.431	0.381	0.127	0.983	1.030	0.224	0.223	NSEI
p_yellow	13	12	2	0.434	0.380	0.132	0.981	1.030	0.226	0.239	NSEI
p_pelagic	13	13	1	0.365	0.357	0.136	0.642	1.057	0.133	0.060	NSEI
p_pelagic	13	13	2	0.158	0.146	0.048	0.346	1.046	0.079	0.153	NSEI
p_yellow	13	13	1	0.433	0.384	0.135	0.980	1.029	0.223	0.219	NSEI
p_yellow	13	13	2	0.432	0.379	0.130	0.985	1.031	0.224	0.236	NSEI
p_pelagic	13	14	1	0.343	0.328	0.139	0.621	1.052	0.124	0.121	NSEI
p_pelagic	13	14	2	0.155	0.141	0.049	0.340	1.076	0.076	0.185	NSEI
p_yellow	13	14	1	0.428	0.370	0.131	0.981	1.028	0.225	0.258	NSEI
p_yellow	13	14	2	0.436	0.384	0.129	0.981	1.027	0.225	0.232	NSEI
p_pelagic	13	15	1	0.315	0.299	0.120	0.578	1.044	0.120	0.133	NSEI
p_pelagic	13	15	2	0.145	0.133	0.046	0.321	1.063	0.072	0.166	NSEI
p_yellow	13	15	1	0.426	0.376	0.124	0.982	1.034	0.224	0.223	NSEI
p_yellow	13	15	2	0.432	0.380	0.131	0.985	1.032	0.224	0.233	NSEI
p_pelagic	13	16	1	0.350	0.341	0.141	0.616	1.038	0.125	0.072	NSEI
p_pelagic	13	16	2	0.162	0.149	0.053	0.350	1.046	0.078	0.166	NSEI
p_yellow	13	16	1	0.429	0.377	0.128	0.979	1.032	0.224	0.233	NSEI
p_yellow	13	16	2	0.431	0.383	0.127	0.983	1.026	0.224	0.215	NSEI
p_pelagic	13	17	1	0.363	0.355	0.159	0.620	1.030	0.123	0.065	NSEI
p_pelagic	13	17	2	0.163	0.151	0.055	0.340	1.062	0.076	0.159	NSEI
p_yellow	13	17	1	0.431	0.381	0.126	0.982	1.033	0.225	0.222	NSEI
p_yellow	13	17	2	0.434	0.382	0.132	0.985	1.029	0.224	0.232	NSEI
p_pelagic	13	18	1	0.368	0.356	0.159	0.631	1.041	0.124	0.097	NSEI
p_pelagic	13	18	2	0.167	0.154	0.056	0.349	1.038	0.076	0.170	NSEI
p_yellow	13	18	1	0.432	0.382	0.121	0.981	1.030	0.224	0.223	NSEI
p_yellow	13	18	2	0.434	0.381	0.132	0.981	1.030	0.224	0.237	NSEI
p_pelagic	13	19	1	0.366	0.353	0.156	0.633	1.036	0.123	0.106	NSEI
p_pelagic	13	19	2	0.170	0.158	0.059	0.355	1.035	0.078	0.153	NSEI
p_yellow	13	19	1	0.428	0.381	0.124	0.980	1.031	0.223	0.211	NSEI
p_yellow	13	19	2	0.434	0.385	0.134	0.982	1.028	0.225	0.218	NSEI
p_pelagic	13	2	1	0.343	0.328	0.129	0.628	1.077	0.133	0.113	NSEI
p_pelagic	13	2	2	0.149	0.132	0.043	0.337	1.083	0.078	0.217	NSEI
p_yellow	13	2	1	0.431	0.382	0.127	0.982	1.033	0.224	0.219	NSEI
p_yellow	13	2	2	0.435	0.386	0.130	0.981	1.028	0.224	0.219	NSEI
p_pelagic	13	20	1	0.380	0.372	0.167	0.651	1.032	0.123	0.065	NSEI
p_pelagic	13	20	2	0.176	0.162	0.061	0.353	1.027	0.076	0.185	NSEI
p_yellow	13	20	1	0.426	0.376	0.124	0.981	1.031	0.224	0.223	NSEI
p_yellow	13	20	2	0.430	0.375	0.126	0.981	1.025	0.225	0.244	NSEI
p_pelagic	13	21	1	0.397	0.390	0.186	0.649	1.031	0.122	0.057	NSEI
p_pelagic	13	21	2	0.185	0.171	0.065	0.377	1.023	0.082	0.170	NSEI
p_yellow	13	21	1	0.428	0.378	0.133	0.983	1.035	0.223	0.225	NSEI
p_yellow	13	21	2	0.431	0.381	0.131	0.984	1.027	0.224	0.223	NSEI
p_pelagic	13	22	2	0.184	0.173	0.064	0.365	1.016	0.078	0.141	NSEI
p_yellow	13	22	1	0.424	0.373	0.122	0.980	1.038	0.224	0.228	NSEI
p_yellow	13	22	2	0.432	0.381	0.128	0.983	1.025	0.223	0.228	NSEI
p_pelagic	13	23	2	0.117	0.108	0.043	0.242	1.016	0.051	0.175	NSEI
p_yellow	13	23	1	0.392	0.335	0.103	0.974	1.030	0.225	0.254	NSEI
p_yellow	13	23	2	0.483	0.446	0.165	0.987	1.018	0.216	0.171	NSEI
p_yellow	13	24	1	0.395	0.339	0.103	0.976	1.035	0.228	0.246	NSEI
p_yellow	13	24	2	0.481	0.441	0.162	0.986	1.021	0.217	0.184	NSEI
p_yellow	13	25	1	0.396	0.341	0.108	0.975	1.031	0.223	0.246	NSEI
p_yellow	13	25	2	0.480	0.441	0.154	0.986	1.024	0.217	0.180	NSEI
p_yellow	13	26	1	0.387	0.330	0.098	0.975	1.027	0.224	0.255	NSEI
p_yellow	13	26	2	0.471	0.433	0.158	0.986	1.027	0.218	0.175	NSEI
p_yellow	13	27	1	0.380	0.326	0.098	0.967	1.026	0.219	0.247	NSEI
p_yellow	13	27	2	0.473	0.436	0.155	0.984	1.023	0.214	0.173	NSEI
p_yellow	13	28	1	0.388	0.342	0.101	0.968	1.025	0.215	0.214	NSEI
p_yellow	13	28	2	0.476	0.440	0.161	0.981	1.029	0.208	0.173	NSEI
p_yellow	13	29	1	0.365	0.318	0.099	0.936	1.024	0.200	0.235	NSEI
p_yellow	13	29	2	0.468	0.442	0.156	0.956	1.020	0.200	0.130	NSEI
p_pelagic	13	3	1	0.344	0.329	0.130	0.634	1.074	0.133	0.113	NSEI
p_pelagic	13	3	2	0.150	0.133	0.042	0.343	1.075	0.079	0.214	NSEI
p_yellow	13	3	1	0.431	0.378	0.127	0.980	1.031	0.224	0.237	NSEI
p_yellow	13	3	2	0.434	0.385	0.130	0.981	1.029	0.225	0.218	NSEI
p_pelagic	13	4	1	0.344	0.334	0.129	0.622	1.054	0.131	0.076	NSEI
p_pelagic	13	4	2	0.149	0.133	0.042	0.349	1.084	0.079	0.203	NSEI
p_yellow	13	4	1	0.433	0.386	0.128	0.983	1.027	0.224	0.209	NSEI
p_yellow	13	4	2	0.439	0.385	0.134	0.984	1.033	0.225	0.240	NSEI
p_pelagic	13	5	1	0.345	0.332	0.125	0.622	1.053	0.131	0.099	NSEI
p_pelagic	13	5	2	0.150	0.134	0.042	0.342	1.054	0.079	0.202	NSEI
p_yellow	13	5	1	0.429	0.378	0.129	0.984	1.032	0.224	0.228	NSEI
p_yellow	13	5	2	0.434	0.380	0.136	0.984	1.033	0.225	0.240	NSEI
p_pelagic	13	6	1	0.350	0.337	0.131	0.627	1.057	0.133	0.098	NSEI
p_pelagic	13	6	2	0.150	0.134	0.042	0.352	1.077	0.079	0.202	NSEI
p_yellow	13	6	1	0.430	0.378	0.128	0.980	1.035	0.226	0.230	NSEI
p_yellow	13	6	2	0.432	0.383	0.131	0.984	1.030	0.225	0.218	NSEI
p_pelagic	13	7	1	0.350	0.340	0.130	0.631	1.053	0.133	0.075	NSEI
p_pelagic	13	7	2	0.153	0.134	0.045	0.345	1.058	0.081	0.236	NSEI
p_yellow	13	7	1	0.431	0.385	0.128	0.982	1.027	0.224	0.206	NSEI
p_yellow	13	7	2	0.438	0.389	0.131	0.982	1.030	0.224	0.218	NSEI
p_pelagic	13	8	1	0.350	0.339	0.131	0.625	1.062	0.132	0.083	NSEI
p_pelagic	13	8	2	0.151	0.135	0.044	0.348	1.072	0.079	0.203	NSEI
p_yellow	13	8	1	0.429	0.379	0.121	0.982	1.031	0.223	0.224	NSEI
p_yellow	13	8	2	0.441	0.395	0.135	0.984	1.033	0.224	0.205	NSEI
p_pelagic	13	9	1	0.349	0.335	0.133	0.631	1.075	0.129	0.108	NSEI
p_pelagic	13	9	2	0.152	0.138	0.045	0.342	1.081	0.077	0.181	NSEI
p_yellow	13	9	1	0.434	0.382	0.130	0.984	1.030	0.224	0.232	NSEI
p_yellow	13	9	2	0.435	0.389	0.127	0.982	1.029	0.224	0.205	NSEI
p_yellow	14	1	1	0.585	0.567	0.256	0.985	1.010	0.189	0.095	NSEO
p_yellow	14	10	1	0.579	0.563	0.254	0.985	1.015	0.189	0.085	NSEO
p_yellow	14	11	1	0.587	0.572	0.254	0.985	1.010	0.191	0.079	NSEO
p_yellow	14	12	1	0.587	0.572	0.266	0.984	1.013	0.189	0.079	NSEO
p_yellow	14	12	2	0.622	0.614	0.287	0.989	1.011	0.185	0.043	NSEO
p_yellow	14	13	1	0.586	0.569	0.251	0.984	1.014	0.188	0.090	NSEO
p_yellow	14	14	1	0.586	0.572	0.258	0.988	1.012	0.190	0.074	NSEO
p_yellow	14	14	2	0.631	0.622	0.303	0.989	1.011	0.180	0.050	NSEO
p_yellow	14	16	2	0.627	0.621	0.299	0.988	1.011	0.182	0.033	NSEO
p_yellow	14	18	1	0.583	0.564	0.252	0.986	1.011	0.192	0.099	NSEO
p_yellow	14	18	2	0.624	0.614	0.285	0.989	1.011	0.184	0.054	NSEO
p_yellow	14	19	1	0.589	0.575	0.252	0.986	1.010	0.191	0.073	NSEO
p_yellow	14	19	2	0.625	0.617	0.288	0.989	1.014	0.181	0.044	NSEO
p_yellow	14	2	1	0.584	0.567	0.262	0.985	1.013	0.188	0.090	NSEO
p_yellow	14	20	1	0.586	0.572	0.251	0.985	1.013	0.190	0.074	NSEO
p_yellow	14	21	1	0.589	0.578	0.255	0.987	1.014	0.189	0.058	NSEO
p_yellow	14	22	2	0.625	0.618	0.296	0.991	1.012	0.181	0.039	NSEO
p_yellow	14	23	1	0.618	0.611	0.299	0.986	1.013	0.178	0.039	NSEO
p_yellow	14	23	2	0.628	0.620	0.284	0.988	1.010	0.183	0.044	NSEO
p_yellow	14	26	1	0.595	0.585	0.264	0.985	1.012	0.188	0.053	NSEO
p_yellow	14	26	2	0.625	0.620	0.287	0.990	1.012	0.184	0.027	NSEO
p_yellow	14	28	1	0.591	0.581	0.272	0.985	1.017	0.182	0.055	NSEO
p_yellow	14	28	2	0.621	0.614	0.284	0.988	1.010	0.184	0.038	NSEO
p_yellow	14	29	1	0.592	0.588	0.262	0.976	1.012	0.180	0.022	NSEO
p_yellow	14	29	2	0.618	0.615	0.282	0.982	1.013	0.180	0.017	NSEO
p_yellow	14	4	1	0.586	0.571	0.254	0.984	1.012	0.190	0.079	NSEO
p_dsr	14	44	1	0.076	0.059	0.015	0.234	1.014	0.058	0.293	NSEO
p_yellow	14	5	2	0.623	0.613	0.274	0.989	1.010	0.185	0.054	NSEO
p_yellow	14	6	1	0.588	0.572	0.264	0.986	1.012	0.189	0.085	NSEO
p_yellow	14	6	2	0.629	0.621	0.306	0.987	1.013	0.180	0.044	NSEO
p_yellow	14	8	1	0.586	0.570	0.254	0.985	1.013	0.189	0.085	NSEO
p_pelagic	16	1	1	0.380	0.368	0.140	0.664	1.062	0.137	0.088	SSEO
p_pelagic	16	1	2	0.359	0.348	0.125	0.629	1.052	0.134	0.082	SSEO
p_pelagic	16	10	1	0.386	0.377	0.157	0.671	1.045	0.135	0.067	SSEO
p_pelagic	16	10	2	0.366	0.357	0.140	0.652	1.040	0.134	0.067	SSEO
p_pelagic	16	11	1	0.386	0.378	0.156	0.663	1.049	0.134	0.060	SSEO
p_pelagic	16	11	2	0.366	0.362	0.140	0.637	1.038	0.133	0.030	SSEO
p_pelagic	16	12	1	0.390	0.379	0.153	0.668	1.044	0.136	0.081	SSEO
p_pelagic	16	12	2	0.364	0.358	0.139	0.636	1.041	0.131	0.046	SSEO
p_pelagic	16	13	1	0.392	0.384	0.156	0.657	1.042	0.134	0.060	SSEO
p_pelagic	16	13	2	0.367	0.358	0.144	0.641	1.039	0.133	0.068	SSEO
p_pelagic	16	14	1	0.375	0.366	0.156	0.651	1.032	0.130	0.069	SSEO
p_pelagic	16	14	2	0.366	0.357	0.142	0.653	1.043	0.132	0.068	SSEO
p_pelagic	16	15	1	0.381	0.374	0.162	0.646	1.029	0.128	0.055	SSEO
p_pelagic	16	15	2	0.369	0.360	0.144	0.635	1.030	0.129	0.070	SSEO
p_pelagic	16	16	1	0.382	0.373	0.153	0.645	1.033	0.131	0.069	SSEO
p_pelagic	16	16	2	0.370	0.360	0.144	0.636	1.042	0.128	0.078	SSEO
p_pelagic	16	17	1	0.427	0.421	0.194	0.688	1.042	0.127	0.047	SSEO
p_pelagic	16	17	2	0.383	0.375	0.154	0.658	1.024	0.129	0.062	SSEO
p_pelagic	16	18	1	0.416	0.410	0.187	0.668	1.028	0.127	0.047	SSEO
p_pelagic	16	18	2	0.386	0.377	0.161	0.648	1.024	0.128	0.070	SSEO
p_pelagic	16	19	1	0.399	0.392	0.185	0.651	1.021	0.123	0.057	SSEO
p_pelagic	16	19	2	0.385	0.379	0.163	0.647	1.027	0.128	0.047	SSEO
p_pelagic	16	2	1	0.381	0.374	0.137	0.664	1.062	0.139	0.050	SSEO
p_pelagic	16	2	2	0.362	0.352	0.132	0.646	1.056	0.136	0.074	SSEO
p_pelagic	16	20	1	0.404	0.396	0.185	0.658	1.016	0.123	0.065	SSEO
p_pelagic	16	20	2	0.388	0.379	0.166	0.648	1.012	0.128	0.070	SSEO
p_pelagic	16	21	1	0.388	0.379	0.172	0.642	1.013	0.123	0.073	SSEO
p_pelagic	16	21	2	0.385	0.379	0.169	0.648	1.013	0.125	0.048	SSEO
p_pelagic	16	3	1	0.379	0.370	0.140	0.662	1.052	0.137	0.066	SSEO
p_pelagic	16	3	2	0.362	0.352	0.130	0.653	1.056	0.136	0.073	SSEO
p_dsr	16	31	1	0.402	0.400	0.375	0.434	1.783	0.016	0.129	SSEO
p_dsr	16	36	1	0.603	0.603	0.573	0.634	1.014	0.015	0.000	SSEO
p_slope	16	36	1	0.060	0.061	0.040	0.078	1.068	0.009	-0.107	SSEO
p_dsr	16	37	1	0.616	0.616	0.587	0.646	1.013	0.015	0.000	SSEO
p_slope	16	37	1	0.065	0.065	0.044	0.083	1.089	0.010	0.000	SSEO
p_pelagic	16	4	1	0.381	0.370	0.144	0.659	1.047	0.135	0.081	SSEO
p_pelagic	16	4	2	0.359	0.354	0.130	0.624	1.051	0.132	0.038	SSEO
p_pelagic	16	42	1	0.833	0.833	0.826	0.840	1.111	0.003	0.000	SSEO
p_yellow	16	42	1	0.407	0.407	0.385	0.428	1.038	0.011	0.000	SSEO
p_dsr	16	42	1	0.451	0.441	0.405	0.509	3.601	0.031	0.326	SSEO
p_slope	16	42	1	0.167	0.180	0.109	0.210	4.888	0.032	-0.405	SSEO
p_pelagic	16	5	1	0.382	0.375	0.148	0.668	1.050	0.138	0.051	SSEO
p_pelagic	16	5	2	0.357	0.349	0.127	0.633	1.059	0.133	0.060	SSEO
p_pelagic	16	6	1	0.385	0.378	0.141	0.658	1.058	0.136	0.052	SSEO
p_pelagic	16	6	2	0.361	0.355	0.133	0.635	1.059	0.133	0.045	SSEO
p_pelagic	16	7	1	0.378	0.375	0.138	0.657	1.047	0.135	0.022	SSEO
p_pelagic	16	7	2	0.362	0.349	0.132	0.637	1.066	0.135	0.096	SSEO
p_pelagic	16	8	1	0.388	0.379	0.146	0.662	1.052	0.135	0.066	SSEO
p_pelagic	16	8	2	0.361	0.354	0.133	0.646	1.057	0.134	0.052	SSEO
p_pelagic	16	9	1	0.384	0.378	0.151	0.657	1.046	0.134	0.045	SSEO
p_pelagic	16	9	2	0.362	0.353	0.136	0.631	1.042	0.133	0.068	SSEO
p_pelagic	4	1	1	0.559	0.564	0.292	0.793	1.029	0.130	-0.039	PWSO
p_pelagic	4	1	2	0.555	0.565	0.270	0.798	1.032	0.134	-0.075	PWSO
p_pelagic	4	10	1	0.565	0.573	0.300	0.794	1.029	0.130	-0.062	PWSO
p_pelagic	4	10	2	0.559	0.565	0.289	0.792	1.030	0.130	-0.046	PWSO
p_pelagic	4	11	1	0.567	0.574	0.309	0.799	1.038	0.129	-0.054	PWSO
p_pelagic	4	11	2	0.557	0.562	0.294	0.793	1.028	0.131	-0.038	PWSO
p_pelagic	4	12	1	0.558	0.558	0.295	0.793	1.030	0.132	0.000	PWSO
p_pelagic	4	12	2	0.554	0.557	0.285	0.795	1.025	0.133	-0.023	PWSO
p_pelagic	4	13	1	0.566	0.574	0.295	0.809	1.034	0.132	-0.060	PWSO
p_pelagic	4	13	2	0.560	0.567	0.295	0.791	1.028	0.130	-0.054	PWSO
p_pelagic	4	14	1	0.561	0.566	0.302	0.795	1.029	0.130	-0.038	PWSO
p_pelagic	4	14	2	0.557	0.562	0.297	0.786	1.023	0.131	-0.038	PWSO
p_pelagic	4	15	1	0.559	0.562	0.292	0.797	1.033	0.131	-0.023	PWSO
p_pelagic	4	15	2	0.561	0.565	0.305	0.804	1.027	0.128	-0.031	PWSO
p_pelagic	4	16	1	0.559	0.564	0.292	0.783	1.021	0.128	-0.039	PWSO
p_pelagic	4	16	2	0.555	0.560	0.289	0.795	1.031	0.130	-0.038	PWSO
p_pelagic	4	17	1	0.562	0.565	0.299	0.800	1.038	0.131	-0.023	PWSO
p_pelagic	4	17	2	0.558	0.564	0.299	0.796	1.027	0.131	-0.046	PWSO
p_pelagic	4	18	1	0.561	0.564	0.300	0.798	1.028	0.130	-0.023	PWSO
p_pelagic	4	18	2	0.558	0.564	0.291	0.800	1.019	0.132	-0.045	PWSO
p_pelagic	4	19	1	0.555	0.558	0.290	0.798	1.027	0.131	-0.023	PWSO
p_pelagic	4	19	2	0.557	0.567	0.300	0.797	1.038	0.132	-0.076	PWSO
p_pelagic	4	2	1	0.562	0.570	0.298	0.792	1.032	0.131	-0.061	PWSO
p_pelagic	4	2	2	0.557	0.562	0.291	0.794	1.030	0.131	-0.038	PWSO
p_pelagic	4	20	2	0.471	0.468	0.244	0.708	1.023	0.121	0.025	PWSO
p_pelagic	4	21	2	0.410	0.409	0.245	0.589	1.017	0.088	0.011	PWSO
p_pelagic	4	22	2	0.465	0.465	0.284	0.652	1.015	0.095	0.000	PWSO
p_pelagic	4	3	1	0.565	0.575	0.285	0.803	1.035	0.132	-0.076	PWSO
p_pelagic	4	3	2	0.555	0.563	0.292	0.794	1.035	0.132	-0.061	PWSO
p_pelagic	4	4	1	0.562	0.568	0.296	0.798	1.017	0.129	-0.046	PWSO
p_pelagic	4	4	2	0.555	0.562	0.287	0.803	1.038	0.134	-0.052	PWSO
p_pelagic	4	5	1	0.562	0.570	0.307	0.792	1.032	0.129	-0.062	PWSO
p_pelagic	4	5	2	0.555	0.559	0.294	0.787	1.022	0.128	-0.031	PWSO
p_pelagic	4	6	1	0.563	0.569	0.307	0.790	1.032	0.128	-0.047	PWSO
p_pelagic	4	6	2	0.560	0.566	0.299	0.795	1.029	0.130	-0.046	PWSO
p_pelagic	4	7	1	0.562	0.571	0.295	0.795	1.022	0.132	-0.068	PWSO
p_pelagic	4	7	2	0.560	0.565	0.294	0.793	1.033	0.130	-0.039	PWSO
p_pelagic	4	8	1	0.561	0.567	0.301	0.795	1.031	0.129	-0.046	PWSO
p_pelagic	4	8	2	0.556	0.561	0.291	0.798	1.038	0.133	-0.037	PWSO
p_pelagic	4	9	1	0.567	0.572	0.297	0.800	1.027	0.129	-0.039	PWSO
p_pelagic	4	9	2	0.560	0.565	0.298	0.801	1.041	0.133	-0.038	PWSO
p_black	8	40	1	0.440	0.412	0.340	0.949	1.126	0.128	0.219	afognak
p_black	8	48	1	0.525	0.491	0.414	0.963	1.053	0.129	0.263	afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable	area_n	year/species	mean	median	lower	upper	Rhat	sd	skew	area
beta0_pH	10	1	0.265	0.292	-0.231	0.616	1.119	0.208	-0.130	northeast
beta1_pH	10	1	2.174	2.058	1.460	4.570	1.258	0.615	0.189	northeast
beta2_pH	10	1	1.706	0.912	0.140	9.578	1.319	2.330	0.341	northeast
beta2_pH	10	3	0.868	0.277	0.105	11.247	1.243	2.335	0.253	northeast
beta3_pH	10	1	32.532	32.273	30.069	38.076	1.182	1.738	0.149	northeast
beta0_pH	11	1	0.103	0.108	-0.551	0.676	1.132	0.310	-0.016	CSEO
beta2_pH	11	1	0.467	0.408	0.221	1.053	1.104	0.226	0.260	CSEO
beta0_pelagic	11		-0.101	-0.093	-0.473	0.196	1.284	0.166	-0.048	CSEO
beta1_pelagic	11		1.088	1.090	0.678	1.531	1.282	0.227	-0.009	CSEO
beta2_pelagic	11		1.870	1.377	0.303	6.032	1.179	1.534	0.321	CSEO
beta0_pelagic	13		-0.728	-0.698	-1.442	-0.019	1.226	0.389	-0.077	NSEI
beta1_pelagic	13		2.368	2.340	1.346	3.395	1.215	0.563	0.050	NSEI
beta3_pelagic	13		31.557	31.037	25.510	39.680	1.184	3.725	0.140	NSEI
beta3_pelagic	14		34.340	34.189	31.378	38.454	1.102	1.739	0.087	NSEO
beta0_pH	15	1	0.086	0.103	-0.333	0.405	1.108	0.186	-0.091	SSEI
beta0_yellow	15		-1.261	-1.235	-1.785	-1.022	1.175	0.174	-0.149	SSEI
beta1_yellow	15		0.175	0.079	0.011	0.906	1.193	0.232	0.413	SSEI
beta2_yellow	15		-5.846	-4.131	-20.009	-0.182	1.330	5.431	-0.316	SSEI
beta0_pH	16	5	0.252	0.240	0.157	0.412	1.136	0.067	0.178	SSEO
beta3_pH	16	4	41.579	41.585	40.888	42.388	1.136	0.402	-0.015	SSEO
beta0_pelagic	16		-0.539	-0.485	-1.399	-0.004	1.230	0.383	-0.141	SSEO
beta1_pelagic	16		1.646	1.510	0.816	3.002	1.175	0.599	0.227	SSEO
beta0_dsr	16		0.160	0.160	0.128	0.193	1.224	0.016	0.000	SSEO
beta0_slope	16		-2.650	-2.651	-2.730	-2.562	1.174	0.043	0.024	SSEO
mu_beta2_pH	2	1	1.762	0.941	0.163	9.219	1.377	2.490	0.330	NG
beta2_yellow	2		-9.401	-9.118	-20.733	-0.522	1.216	5.790	-0.049	NG
beta1_pH	4	1	1.960	1.882	1.246	3.700	1.261	0.535	0.146	PWSO
beta2_pH	4	1	1.029	0.790	0.194	2.896	1.185	0.741	0.323	PWSO
beta0_pelagic	4		0.268	0.293	-0.292	0.702	1.214	0.243	-0.103	PWSO
beta1_pelagic	4		0.778	0.755	0.267	1.379	1.202	0.274	0.084	PWSO
beta0_yellow	4		0.844	0.949	-0.124	1.188	1.138	0.333	-0.315	PWSO
beta2_pH	5	1	2.357	1.570	0.282	10.127	1.346	2.662	0.296	BSAI
beta0_pH	6	1	-0.449	-0.329	-1.603	0.030	1.253	0.405	-0.296	SOKO2SAP
beta1_pH	6	1	1.933	1.788	1.080	3.367	1.101	0.621	0.234	SOKO2SAP
beta2_pH	6	1	1.644	0.577	0.141	10.824	1.367	2.816	0.379	SOKO2SAP
beta3_pH	6	1	33.025	33.148	25.836	38.926	1.137	2.871	-0.043	SOKO2SAP
beta1_pH	8	1	2.085	1.944	1.076	3.838	1.151	0.703	0.200	afognak
beta2_pH	8	1	1.363	0.465	0.150	9.360	1.436	2.309	0.389	afognak
beta2_pH	8	3	1.976	1.135	0.143	11.512	1.174	2.524	0.333	afognak
beta1_pH	9	1	1.963	1.887	1.240	3.078	1.147	0.471	0.161	eastside
beta2_pH	9	1	1.437	0.621	0.188	8.572	1.352	2.235	0.365	eastside
beta2_pH	9	3	0.761	0.213	0.092	11.002	1.269	2.335	0.235	eastside
lambda_H	5		2.908	0.621	0.026	21.835	1.139	8.713	0.262	BSAI

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic	non.BSAI.SOKO.Harv	BSAI.SOKO.Harv	non.BSAI.SOKO.Rel	BSAI.SOKO.Rel
n unconverged at Rhat = 1.01	811.00	274.00	1,215.00	347.00
Prop. of total	0.75	0.25	0.78	0.22
Prop. median less than 1	0.17	0.19	0.28	0.41
Prop. median less than 10	0.24	0.47	0.44	0.70
Prop. median less than 50	0.34	0.69	0.60	0.78
Median of the median values	402.27	12.24	18.77	1.90
Mean of the median values	1,437.40	561.68	452.16	187.14
Mean skew of posteriors	0.20	0.28	0.26	0.26

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic	non.BSAI.SOKO.Harv	BSAI.SOKO.Harv	non.BSAI.SOKO.Rel	BSAI.SOKO.Rel
n unconverged at Rhat = 1.1	12.00	11.00	59.00	38.00
Prop. of total	0.52	0.48	0.61	0.39
Prop. median less than 1	0.25	0.55	0.64	0.66
Prop. median less than 10	0.25	1.00	0.69	0.84
Prop. median less than 50	0.25	1.00	0.83	0.89
Median of the median values	286.28	0.98	0.01	0.09
Mean of the median values	635.64	1.12	87.35	19.43
Mean skew of posteriors	0.22	0.21	0.15	0.18

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable	Freq
beta0_dsr	1
beta0_pelagic	4
beta0_pH	4
beta0_slope	1
beta0_yellow	2
beta1_pelagic	4
beta1_pH	4
beta1_yellow	1
beta2_pelagic	1
beta2_pH	9
beta2_yellow	2
beta3_pelagic	1
beta3_pH	3
H_ay	3
H_ayg	1
H_ayu	3
Hb_ay	3
Hb_ayg	2
Hb_ayu	3
Hd_ay	1
Hd_ayg	2
Ho_ay	1
Ho_ayg	1
Ho_ayu	3
Hp_ay	3
Hp_ayg	1
Hp_ayu	3
Hs_ay	2
Hs_ayg	2
Hs_ayu	1
Hy_ay	4
Hy_ayg	2
Hy_ayu	5
lambda_H	1
mu_beta2_pH	1
p_black	1
p_dsr	2
p_pelagic	1
p_slope	2
pH	16
R_ay	3
R_ayg	3
R_ayu	5
Rb_ay	3
Rb_ayg	4
Rb_ayu	3
Rd_ay	5
Rd_ayg	6
Rd_ayu	2
re_dsr	2
re_pelagic	23
re_pH	22
re_slope	4
Ro_ay	19
Ro_ayg	21
Ro_ayu	32
Rp_ay	3
Rp_ayg	3
Rp_ayu	3
Rs_ay	5
Rs_ayg	11
Rs_ayu	16
Ry_ay	6
Ry_ayg	5
Ry_ayu	13
sd_a_wt	13
sd_wt	14

Parameter estimates:

Summary Table of Parameter Estimates
Parameter	mean	sd	Lower_CI	Median	Upper_CI

Rockfish Historical Harvest Estimates

Phil Joy & Adam Reimer

February 2025